Band structure and optical properties of opal photonic crystals 
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A theoretical approach for the interpretation of reflectance spectra of opal photonic crystals with 
fee structure and (111) surface orientation is presented. It is based on the calculation of photonic 
bands and density of states corresponding to a specified angle of incidence in air. The results yield 
a clear distinction between diffraction in the direction of light propagation by (111) family planes 
(leading to the formation of a stop band) and diffraction in other directions by higher-order planes 
(corresponding to the excitation of photonic modes in the crystal). Reflectance measurements on 
artificial opals made of self- assembled polystyrene spheres are analyzed according to the theoretical 
scheme and give evidence of diffraction by higher-order crystalline planes in the photonic structure. 

PACS numbers: 42.70.Qs, 41.20.Jb, 78.40.-q 



I. INTRODUCTION 

Three-dimensional photonic crystals are intensively 
studied with the goal of achieving control of light propa- 
gation and, especially, a complete photonic band gap in 
all directionsi2i2i4 Among the possible structures, fec- 
based systems like opals and inverse opals are of great 
interest as they can be produced with bottom-up ap- 
proaches based on casting from colloidal solutions, possi- 
bly followed by infiltration and template removali^Si^ 
The direct opal structure (e.g., polystyrene or silica 
spheres in air) has weak effects on the total photonic 
density of states (DOS) but does produce a well-defined 
stop band in the direction of light propagation. If the 
structure crystallizes in an fee lattice with (111) surface 
orientation, as is usually the case, normal-incidence re- 
flectance or transmittance probe the photonic gap in the 
TL direction of the fee Brillouin zone. Instead, the inverse 
opal structure (e.g., air spheres in TiC>2 or in silicon) pro- 
duces a strong effect on the photonic DOS and may even 
result in a complete photonic band gap if the dielectric 
constant and the filling fraction of the infiltrated material 
are large enoughAiS 

The optical properties of fee colloidal pho- 
tonic crystals have been the subject of numerous 
c-pcrimcnta l 11 ' 12 - 13 - 14,15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24,25 - 26 - 27 - 28 - 29 
and theoretical 9 i 10 i 20 i 34 i 35 i 36 i 37 i 38 i 39 studies. Many ex- 
perimental investigations are concerned with properties 
of the lowest-order gap, in particular the attenuation 
studied by normal-incidence transmission 1 ^^ and the 
evolution of the stop band in reflectance/transmittance 
spectra at varying angles of incidencei 4 * 1 ^* 1 ^* 2 ^ 2 ^^^ 
and as a function of polarization: 3 *" 1 Detailed studies of 
the high-energy region at normal incidence have recently 
appeared^ The role of disorder in Bragg scattering and 
gap formation has been addressed with measurements 
on multiple-domain^^ 9 " and single-domain crystalsi^Si 2 ^ 



Bragg diffraction has also been studied in a transmission 
geometry by Kossel line analysis^ The change in 
effective refractive index close to the first-order gap and 
in connection with flat bands at high energies has also 
been studied pi&i^i Laser dyes in opal photonic crystals 
show effects related to modified photonic density of 
states in emission spectra 19,23 '— and to interaction of 
photonic modes with excitonsi 4 ^ On the theoretical side, 
the formation of a complete band gap in the inverse 
opal structur o 9 ! 10 and its fragility against disorder 35 
are now well understood. The interplay between 
photonic bands and reflection/transmission properties 
has been addressed both with calculations of optical 



spectr; 



,20,29.34 



and superprism effects^ as well as with 



symmetry analysis^ More recent works focus on the 
effects of stacking faults on the optical properties. 3 ''' 39 

A question of particular interest is the role of diffrac- 
tion from lattice planes other than the (111) family in 
connection with the optical properties at energies higher 
than the first stop band along the TL direction. The 
analysis of Kossel lines in transmission has led to a deter- 
mination of the photonic gap positions for the two lowest 
band gapsi^: the results for the band edges are reported 
as a function of incidence angle along the LW direction 
of the fee Brillouin zone. This interpretation has been 
30,3i.32.,» onC( j m R e f.l34lon the basis of calculations of trans- 



mission spectra and projected photonic bands for given 
values of the wavevector k|| parallel to the crystal sur- 
face. Angle-resolved reflectance measurements give evi- 
dence of multiple diffraction peaks with small splittings 
at large values of the angle of incidence ^iS 5 ^^ which 
have been interpreted as arising from multiple Bragg- 
wave coupling 24 ' 25 - 32 or photonic band branching, 2( i both 
concepts being related to an anticrossing of photonic lev- 
els close to the U or W points of the Brillouin zone. Re- 
lated observations but a different interpretation is given 
in Refl3lL in which additional peaks in reflectance (or 
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FIG. 1: Energy bands for an opal photonic crystal (fee lat- 
tice of close-packed dielectric spheres with ei = 2.53 in air): 

(a) along symmetry directions in the whole Brillouin zone, 

(b) along symmetry directions on the hexagonal face of the 
Brillouin zone. 




FIG. 2: Symmetry points in the Brillouin zone for a fee lattice. 



dips in transmittance) at near-normal incidence are at- 
tributed to a large effective refractive index arising from 
flat bands. The main issue in the interpretation of all 
these results is that spectral features in optical measure- 
ments cannot be unambiguously associated with the pres- 
ence of a photonic gap and/or with diffraction. For re- 
flectance measurements, in particular, a spectral feature 
may follow from the existence of a photonic gap (as is the 
case for the lowest stop band at near-normal incidence) 
but also from the excitation of a photonic mode in the 
crystal: establishing a clear distinction between the two 
effects is the main purpose of the present paper. 

In this work we present a theoretical approach to the 
optical properties of 3D photonic crystals which yields 
further insight in the phenomenon of diffraction from dif- 
ferent families of lattice planes in relation to the presence 
of photonic gaps or photonic bands. We calculate the 
photonic bands and density of states corresponding to a 
fixed external angle of incidence: the DOS is decreased 
when a photonic gap along the propagation direction oc- 
curs, while it increases and shows resonant structures 
when diffraction in directions other than that of specu- 
lar reflection takes place. Calculations of reflection and 
diffraction spectra show that peaks in the DOS corre- 
spond to the onset of diffraction, as well as to additional 
reflectance structures at large values of the angle of inci- 
dence. Reflectance measurements on opal photonic crys- 
tals are presented and interpreted with the above con- 
cepts. In particular, weak features at large angles of in- 
cidence give experimental evidence for the excitation of 
photonic modes in the crystal, i.e., for diffraction from 
planes other than those of the (111) family. 

The results of this work are related to the surface 
reflectance technique for measuring the photonic band 



dispersion. This approach was first applied to photonic 
crystal waveguides, where it led to a determination of the 
dispersion of photonic modes above the light line via the 
observation of resonant structure in reflectance and of 
their dispersion with incidence angle 41424141 The same 
technique of variable-angle reflectance from the crystal 
surface can be applied to the case of macroporous silicon, 
which has a 2D photonic structure and can be considered 
as homogeneous in the direction of the pores, and in this 
case it yields the dispersion of photonic bands in the 2D 
planed The present work shows that, in the case of pho- 
tonic crystals with refractive index modulation in 3D, 
Bragg diffraction off the direction of light propagation 
also produces spectral structures in reflectance, whose 
angular evolution is related to the photonic dispersion 
in a plane parallel to the crystal surface. Preliminary 
results have been recently presented4& 

The rest of this paper is organized as follows. In Sec. II 
we present the theoretical approach, namely the calcula- 
tion of photonic DOS and photonic bands corresponding 
to a specified angle of incidence. In Sec. Ill we describe re- 
flectance experiments on opal photonic crystals and com- 
pare the theoretical predictions with the optical spectra. 
In Section IV the findings of this work are summarized. 



II. PHOTONIC BANDS AND DENSITY OF 
STATES 



We calculate the photonic bands starting from the 
second-order wave equation for the magnetic field 



V x 



1 



V x H 



c z 



(1) 



d 

I 

O 
Q 






0.5 0.6 0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 1 
coa/(27i c) coa/(27i c) coa/(27i c) 

FIG. 3: Reduced density of states for several values of the incidence angle from 8 — 5° (bottom curve) to 9 — 80° (top curve) 
in steps of 5°. The three figures correspond to different values of the azimuthal angle (j>: (a) cj> = 0° + 2rmr/6 (m e 2), or LW 
direction; (b) (j> = 30° + 2mn/3, or LU direction; (c) cj> = —30° + 2mn/3, or LK direction. 



We adopt the plane-wave expansion technique to trans- 
form Eq. (1) into an eigenvalue problem, and calculate 
the inverse dielectric constant by using the procedure of 
Ho, Chan and Soukoulis4£ In Fig. 1 we show the pho- 
tonic bands for an opal with dielectric constants ei = 2.53 
(spheres) 48 i 49 and = 1 (voids), i.e., polystyrene spheres 
in air. The fee Brillouin zone with high symmetry points 
is shown in Fig. 2. The bands are shown in dimensionless 
units uja/(2irc), where a is the lattice constant, along the 
symmetry lines of the full fee Brillouin zone (Fig. la) as 
well as along symmetry lines on the (111) hexagonal face 
(Fig. lb): the latter will be useful for comparing with 
density of states and reflectance spectra at a finite angle 
of incidence. 

The photonic bands in Fig. 1 do not exhibit any com- 
plete gap in the full Brillouin zone (BZ). However, par- 
tial gaps do appear along some symmetry lines of the 
BZ; among these, the lowest energy gap opens along TL 
around oja/(2irc) ~ 0.6. Partial gaps behave as stop 
bands in reflectance spectra when they occur in the di- 
rection of propagation. Photonic modes along the propa- 
gation direction also contribute to the optical properties 
and may give rise to additional features in reflectance, as 
discussed in the rest of this Section. 

In order to distinguish reflectance maxima due to stop 
bands from features which can instead be associated to 
photonic modes, we need first of all to calculate the pho- 
tonic band structure and the density of states (DOS) 



along the direction of propagation, i.e., for a specified 
value of the angle of incidence in air. In the following we 
will refer to these bands and DOS as reduced bands and 
reduced DOS, respectively. The surface of a fee colloidal 
crystal is known to be a (111) crystallographic plane, 
therefore angle-resolved reflectivity experiments are usu- 
ally performed as a function of the angle of incidence with 
respect to a (111) surface. We adopt the same conven- 
tion in our calculations and define a reference frame as 
follows: z is the [1,1,1] or TL direction, x is the [1,0,-1] (a 
LW) direction, and y is the [1,-2,1] (a LK) direction. We 
can therefore express the momentum of an incident pho- 
ton as follows: k l = ^(sin#cos</>,sin#sin</>, cos#), where 
?uv is the energy of the incident photon. With our choice 
of axes, 9 — 0° corresponds to normal incidence (rL di- 
rection) and the azimuthal angle (f> is measured from the 
LW direction. We stress that 9 is the angle of incidence in 
air (which differs from the angle of propagation inside the 
opal, 37 the latter being dependent on the effective refrac- 
tive index and on the presence of additional propagating 
modes when diffraction inside the material occurs). 



The reduced bands and state densities are obtained 
from the full band structure (Fig. 1) by using conserva- 
tions laws: (i) parallel momentum conservation modulo 
a reciprocal lattice vector, and (ii) energy conservation. 
For the photonic bands, these two conditions may be 
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written as 

k|| + G|| = — (sin 9 cos 0, sin 9 sin</>, 0), (2) 
e n (k||,k z ) = hu>, (3) 

where k|| = k| is the component of k parallel to the (111) 
surface, G|| is the parallel component of any reciprocal 
lattice vector G, k z is the component of k parallel to 
TL, and e n (k||,k z ) is the n-th energy band of a photon 
with momentum k. We point out that k|| is fixed by 
the experimental configuration. On the contrary k z and 
the reduced bands are obtained from the solutions of the 
system given by Eqs. J2J,©, which are found by a self- 
consistent procedure. The reduced DOS is obtained from 
the reduced bands as follows: 

N(hu)= ]T ^,(k|,k z )-H (4) 

n,k->0 

where ky is defined by Eq. (0. The above N(hui) is a 
one-dimensional density of states, because it is summed 
only over k z iSi The DOS obtained in this way (convo- 
luted with a Lorentzian of width comparable to a typical 
experimental resolution Au> = 10~ 2 • 2irc/a and averaged 
over photon polarizations) is shown in Fig. 3 for several 
values of 9 and some values of (f> (directions LW, LU and 
LK)£i 

Let us analyze first the results for 9 = 0°. In this case 
the solution of Eqs. J2J and © can be found very simply 
because k|| = 0, and it yields the photonic bands along 
the TL direction. Thus we can directly compare the 9 = 
0° DOS in Fig. 3 with the bands along TL in Fig. la. The 
lowest energy gap shown in Fig. 1 around uja/{2irc) = 0.6 
appears clearly in the reduced DOS at the same values 
of the energy. It arises from diffraction in the backwards 
direction by planes of the (111) family. The photonic 
bands at energies above the gap in Fig. la correspond to 
a photon momentum k outside the first Brillouin zone; 
the bands are folded into the first BZ by subtracting from 
k a reciprocal lattice vector (27r/et)(l, 1,1). These folded 
bands give rise to a finite DOS above the band gap in 
Fig. 3. 

The reduced DOS reflects, of course, the energy dis- 
persion of the photonic bands along the propagation di- 
rection. The latter is almost linear at low frequencies, 
and thus gives rise to a nearly constant one-dimensional 
DOS. For the case shown in Fig. la, the bands along TL 
just above the gap at wa/(2-7rc) = 0.6 exhibit also an al- 
most linear dispersion, thus the corresponding DOS in 
Fig. 3 remains flat above the gap. However, in general, 
the folded bands may exhibit a non-linear dispersion and 
can also show critical points (maxima, minima, or sad- 
dle points). In three dimensions, a critical point results 
in a divergence of the slope of the DOS (van Hove sin- 
gularity). In the present one-dimensional case the DOS 
itself diverges at a van Hove singularity. The very flat 
bands at energy around Loa/{2i:c) = 1.1 (the fifth and 
sixth bands of the photonic structure) show one of such 



critical points at L; as a consequence, the corresponding 
reduced DOS at normal incidence exhibits a Van Hove 
singularity around wa/(2TTc) = 1.1. This critical point is 
out of the scale of Fig. 3 (it is shown in Fig. 4 below); 
however, as we discuss in next paragraph, it moves to 
lower energies when moving from normal incidence. 

The complex, nonlinear dispersion of the photonic 
bands results from the mixing of plane waves with differ- 
ent reciprocal lattice vectors (in directions other than the 
(111) one) in the photonic Bloch mode, thus it implies 
the occurrence of diffraction in directions other than that 
of beam propagation. Notice that this phenomenon can 
occur at any angle of incidence, even for = 0°, when 
the conservation of parallel crystal momentum expressed 
by Eq (2) can be satisfied with G|| =0. For energies be- 
low Lua/(2TTc) = 1.5, the photonic bands are determined 
by reciprocal lattice vectors of lower modulus, the rel- 
evant ones being (000), the eight equivalent {111} and 
the six equivalent {200} ones. Thus the directions of 
diffraction out of the propagation direction (forward or 
backward) are determined by planes of equivalent {111} 
families that are not parallel to the (111) crystal surface, 
as well as by {200} families. At near-normal incidence, 
diffraction directions just above the diffraction cutoff are 
determined by (111), (111), (HI) family planes: on the 
other hand, {200} families become important at large an- 
gles of incidence, as discussed below. We conclude that 
the presence of a van Hove singularity in the reduced 
DOS marks the onset of diffraction from lattice planes 
other than those of the (111) family. 

Let us analyze now the case of finite incidence angles. 
Figure 3 shows how the stop band depend strongly on 
the incidence angle, by moving to higher energies with 
increasing 9. On the other hand, it depends only weakly 
on the azimuthal angle (j). Figure 3 shows also the peaks 
associated to the Van Hove singularity. The latter is at 
uja/(2irc) = 1.1 at normal incidence, but it moves at 
lower energies as 9 increases. This is shown clearly in 
Figs. 3a and 3c, in which this singularity appears as a 
strong peak in the DOS: this peak merges with the gap 
around 9 = 70° and 60°, respectively. Notice that since 
both the gap and the van Hove singularities occur on the 
hexagonal face of the fee Brillouin zone, their angular 
dependence is related to the wavevector dependence of 
the bands shown in Fig. lb. In particular, along the LW 
orientation, the merging of the van Hove singularity with 
the gap occurs close to the W point (see Fig. lb), where 
multiple Bragg scattering from {111} and {200} family 
planes occurs. We point out that the position of the Van 
Hove singularity depends strongly on 0, as we may see 
by comparing Figs. 3a,c with Fig. 3b. In particular, for 
<f> = 30° (LU direction) the Van Hove singularity remains 
almost outside the energy window. 

Photonic modes give rise to several other features in 
the reduced DOS. All these features are associated with 
folded bands and exhibit a similar behavior as a function 
of 9. This can be seen clearly in Fig. 3b. In this case the 
van Hove singularity remains very high in energy for any 
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FIG. 4: (a) Reduced bands, (b) reduced DOS, (c) reflectance and (d) diffraction spectra for <f> = 0° (LW direction) and three 
values 9 — 5°, 40° 70° of the angle of incidence in air. The bands in (a) are plotted in the interval (0, k m ), where k m = VSiv/a 
is the T-L distance. In the limit 9 = the point k m is identical to L. The bands are displayed only for k z inside the first 
Brillouin zone. 



value of 9. The steps seen in the DOS are instead associ- 



ated with bands which have an almost linear dispersion; 
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a step indicates that a new band is becoming accessible. 
Furthermore, Fig. 3b shows also a second photonic gap 
at higher energies, which is present but much weaker for 
4> = 0° (LW direction, Fig. 3a), and is totally absent for 
cf> = -30° (LK direction, Fig. 3c). 

In order to better clarify the ^-dependence of the stop 
bands and of the van Hove singularities, in Fig. 4a, b we 
show the reduced bands and DOS in more detail. The 
bands are shown as a function of k z and for three values 
of 8. We focus on the LW direction (<f> = 0°), but exami- 
nation of the results for other orientations leads to similar 
conclusions. We notice that, since k z is defined only in- 
side the first BZ, its maximum length becomes shorter 
than ja for large 9 and large frequencies (Fig. 4a for 
9 = 70°), due to the truncated shape of the BZ in Fig. 2: 
this occurs when k|| + G|| is outside the hexagonal face 
containing the L point. 

For near-normal incidence, the bands along the direc- 
tion of propagation (9 = 5°) are simply the photonic 
bands along TL, with the gap around u)a/(2irc) ~ 0.6 
and the flat bands around tua/(2'Kc) ~ 1.1. The panels 
in Fig. 4a show that, when 9 increases, the gap moves 
to higher energy, while the flat bands and the associ- 
ated van Hove singularity move downwards to lower en- 
ergy. Eventually the gap and the flat bands overlap for 
9 ~ 70° and ua/(2irc) ~ 0.75, corresponding to the W 
point of the fee Brillouin zone (see Fig. lb) . These fea- 
tures are precisely reflected in the DOS, as we may see in 
the panels of Fig. 4b. In Fig. 4c, d we show reflectance and 
diffraction spectra for the same values of the angle of in- 
cidence, calculated by the scattering-matrix methodi^S 
Diffraction (D) spectra are obtained as D = 1 — R — T 
in terms of reflectance (R) and transmittance (T). The 
stop band appears clearly in reflectance spectra. The on- 
set of diffraction is around wa/(27rc) = 1.1 for 9 = 5°, 
tva/(2nc) = 0.82 for 9 = 40° and wa/(27rc) = 0.75 for 
9 = 70°. The first peak in diffraction spectra correlates 
very clearly with the peak in the reduced DOS (Fig. 4b) 
arising in turn from the flat bands and the Van Hove 
singularity (Fig. 4a). 

We notice that other diffraction peaks and van Hove 
singularities are present in Fig. 4b, d for frequencies above 
the diffraction cutoff. We do not pursue an anal- 
ysis of these additional features, since they have no 
well-resolved counterpart in either theoretical reflectance 
spectra (Fig. 4c) or in the experimental ones (see next Sec- 
tion), and also because the approximation made in the 
theoretical calculation becomes less accurate on increas- 
ing the energy^ We emphasize that diffraction features 
become evident in reflectance spectra at high values of 
the angle of incidence, as clearly seen for 9 — 70° in 
Fig. 4c and discussed below in connection with the exper- 
imental data. 

Thus the analysis of the reduced bands and DOS cor- 
responding to a specified value of the angle of incidence 
leads to a clear distinction between diffraction effects 
(high DOS regions) from gaps or pseudogaps (zero or 
low DOS). In addition, critical points may give rise to van 



Hove singularities which strongly enhance the DOS in the 
associated spectral region. The presence of a van Hove 
singularity marks the excitation of a photonic mode in 
the crystal and corresponds to the occurrence of diffrac- 
tion in directions other than that of propagation. Van 
Hove singularities appear clearly in diffraction spectra 
and become evident also in reflection spectra at large 
values of the angle of incidence, i.e. close to the W 
point, where Bragg scattering is determined by planes 
of both {111} and {200} families. These conclusions 
are important for interpreting reflection experiments on 
polystyrene opals, as discussed in the next Section. 



III. OPTICAL SPECTRA 

A. Experimental 

Bulk polystyrene opals have been grown by monodis- 
perse sphere suspensions in water (Duke Scientific, di- 
ameter 222, 260, 300, 340, 426 nm, refractive index 
1.59)42^ Small vessels with area of about one squared 
centimeter are obtained by sealing on a glass window 
polystyrene or silicone tubes. These vessels are filled with 
the sphere suspension and placed in a wet environment 
to reduce the evaporation rate. The opals have a menis- 
cus shape and a glossy surface showing bright colored 
reflections. 

Scanning Electron Microscopy (SEM) of a 222 nm opal 
surface has been performed with a Leo Stereoscan 440 
(LEO Electron Microscopy Ltd) electronic microscope. 
The surface is very flat on a large scale (tens of microns) 
and shows a regular arrangement of the spheres in a tri- 
angular lattice as expected for their three-dimensional 
packing in a fee lattice when observed along the (111) 
direction. Figure 5 shows a SEM image of 18x12 (^,m) 2 
opal surface. This image can be considered as repre- 
sentative of the most disordered part of the opal surface 
investigated by the optical measurements. This morphol- 
ogy is very different from that of opal films we previ- 
ously investigated, which were characterized by ordered 
domains of spheres of a few square microns areas with 
both triangular and square lattices, tilted with respect 
to each other and separated by disordered zones or deep 
dislocationsi£i£ii For the present bulk opal samples, no 
square lattices are observed and the surface appears as a 
large triangular lattice of spheres with only a few defects. 

Unpolarized variable-angle specular reflectance has 
been measured in the 0.4-4 eV spectral range by means 
of a Fourier- Transform spectrophotometer Bruker IFS66 
with 8 cm -1 spectral resolution. The light of a Xe arc- 
lamp was collimated and then focused to a spot of 100 jim 
diameter on the sample surface. The sample was placed 
on a home-made 9 — 29 goniometer, that allows the inci- 
dence angle to be varied between 5 and 80 degrees, with 
an angular resolution of 1° determined by the angular 
divergence of the incident light cone. Sample visualiza- 
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FIG. 5: SEM Image of the opal surface (18x22 /im 2 



tion and precise allocation of the light spot on the surface 
was achieved by an optical microscope mounted on the 
goniometer. An Ag mirror was used as a reference. The 
overall error bar of the experimental set-up on the abso- 
lute value of the reflectance is within 2%. 

Particular care has been taken in the alignment of the 
sample. Since the top quality surface of the investigated 
sample is located close to one of the sharp border, we 
used the latter as a macroscopic reference to identify 
the crystallographic axes of the surface. As a matter of 
fact, atomic force microscope images recorded at differ- 
ent magnifications showed that the LW crystallographic 
direction has a well defined offset angle with respect to 
the sharp border. Therefore the sample was mounted on 
the microrcflectometcr by aligning the incidence plane 
with the sharp border and than rotating it by the offset 
angle. 



B. Results and discussion 

We apply our theoretical method to understand the 
reflectance spectra for the sample shown in Fig. 5. We 
therefore calculate the reduced bands for a close packed 
fee crystal with sphere diameter of 222 nm (fee lattice 
constant a = 314 nm). We then obtain the reduced DOS 
and calculate the energy position of the low energy gap 
and its dependence on 8, as well as the two most im- 
portant high energy features, i.e., two peaks that are 
associated with photonic modes. The energy position 
of the gap is obtained directly from the reduced bands 
(see Fig. 4), while the position of the peaks is obtained 
from the DOS shown in Fig. 3. The results are shown in 
Fig. 6a and 6c for the directions LW and LK, respectively. 
The angular dependence of these energy positions corre- 
lates with the dispersion of the photonic bands along the 
hexagonal surface of the Brillouin zone (Fig. lb), which 
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FIG. 6: (ac): Theoretical positions of the gap edges (squares) 
and of the lowest energy peak (stars) for the LW (a) and the 
LK (c) direction. The results are for a close-packed fee struc- 
ture of polystyrene spheres (e = 2.53) with a lattice constant 
a = 314 nm. (bd): Experimental reflectance data for sev- 
eral values of 6 = jA9, with A9 = 5° for the LW (b) and 
the LK (d) direction; curves from left to right correspond to 
j — 1, ... ,16 and are shifted by the value 0.2(j — 1). 



gives the energies of Van Hove singularities. Notice that 
the Van Hove singularity for the LK direction has a lower 
energy (at a given value of the incidence angle) as com- 
pared to the LW direction. The Van Hove singularity 
crosses the band gap at 6 = 70° for LW and 6 = 60° 
for LK. The energies of these structures may be directly 
compared with the features observed in the reflectance 
spectra (Fig. 6b and Fig. 6d). 

Before starting a comparison of the theoretical and 
experimental results, we give a description of the lat- 
ter. Figure 6b and 6d show the variable- angle micro re- 
flectance spectra of the opal for 9 ranging from 5° to 80° 
in steps of 5° and along LW and LK, respectively. These 
spectra arc offset horizontally for a more clear view. At 
near-normal incidence a single, intense (reflectance al- 
most identical to that of the Ag reference mirror) , asym- 
metrical band associated with the stop band is observed 
at 2.47 cV with a FWHM of about 0.15 eV. 

When the incidence angle is increased from normal in- 
cidence up to 45°, the reflectance maximum shifts to- 
wards higher energies slightly reducing its intensity and 
modifying its shape. This behavior is accounted for by 
the theoretical band edge positions summarized by the 
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FIG. 7: Enlarged view of experimental reflectance spectra 
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squares in Fig. 6a (LW) and Fig. 6c (LK). The agree- 
ment between calculated and measured gap position as 
well as its dependence on 8 is very satisfactory. 

A more complex situation occurs for higher values of 
the incidence angle. 

We notice that for 6 > 50° the shape of the spectrum 
changes in a substantial way and additional features oc- 
cur. The details of these features depend on the observed 
direction. To give a better idea of the complicated struc- 
ture of the spectra recorded for large 9 values and in order 
to find evidence of the photonic modes, we enlarge the 
energy scale and show the relevant experimental curves 
in Fig. 7a (LW) and Fig. 7b (LK). A structured band 
is observed to regularly shift its maximum in the range 
from 3 to 3.15 eV retaining the dispersion of the stop 
band. It may roughly be assigned to the interplay of the 
stop band and of the excitation of the photonic modes 
observed in the DOS (Fig. 6a and Fig. 6c), which takes 
place for large 9. However we also notice the appearence 
of weak features in the low-energy side of these bands. 
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units, as a function of sphere diameter: stop band central po- 
sition (triangles) and stop band edges (bars) at near-normal 
incidence; Van Hove singularities at 9 = 70° for LK (circles) 
and LW (squares) orientations. 



In particular, for the LW direction the weak feature is at 
about 2.9 eV, while for LK it is at about 2.7 eV. When 
9 increases, these structures shift towards lower energies, 
i.e., they behave like the calculated van Hove singular- 
ity. Spectral features showing similar dispersions have 
been previously observed in the transmission spectra of 
polystyrene^* and silica^* opal films. The low energy 
calculated peak displayed in Fig. 6a and Fig. 6c for high 
incidence angles is also reported as an arrow in Fig. 7a 
and Fig. 7b to show a direct comparison with the ex- 
perimental spectra for each value of 9. We find a very 
good agreement between the calculated peaks and the 
weak structures observed in the experimental spectra for 
both LW and LK directions. Moreover, the possibility of 
observing a signature of diffraction effects in reflectance 
spectra has been proved by the calculated spectra shown 
in Fig. 4c. We therefore assign the features having disper- 
sion opposite to the stop band to van Hove singularities 
in the reduced DOS, i.e., to diffraction from family planes 
other than those of the (111) one. For large values of the 
incidence angle along the LW orientation, multiple Bragg 
scattering close to the W point takes place and diffrac- 
tion inside the opal is determined by {111} as well as by 
{200} family planes. 

In order to further support the assignment of the spec- 
tral features observed in R spectra with van Hove singu- 
larities, we checked their scaling with the lattice constant 
(related to the sphere diameter). Figure 8 shows the ex- 
perimental spectral position for the stop band maximum 
and its width at near normal incidence, as well as van 
Hove singularities measured along the LW and LK direc- 
tions at 9 = 70° for opals with different sphere diameters. 
All these energies, expressed in dimensionless units, are 
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clearly independent of the lattice constant thus unam- 
biguously showing that the spectral structures are due 
to intrinsic properties of the photonic crystal band struc- 
ture. 

It should be remarked that although the presence of 
fiat bands is essential in giving rise to van Hove singu- 
larities of the DOS and to diffraction features, the inter- 
pretation of the present experiments goes beyond that 
of Refl3ll in which related observations are attributed 
only to a higher effective refractive index arising from 
flat bands and diffraction is not considered. Indeed, when 
diffraction starts to occur, the optical properties of a pho- 
tonic crystal cannot be described by a single effective in- 
dex as several modes propagate in the crystal at the same 
frequency. The present interpretation does not contra- 
dict the multiple Bragg-wave coupling or band branching 
mechanisms previously introduced^^iSSiS and it con- 
firms the identification of multiple diffraction as arising 
from {111} or {200} family planes for the different orien- 
tations. Extending the previous works, our observations 
indicate that diffraction-related features in reflectance 
spectra become increasingly visible for higher values of 
the angle of incidence: i.e., the onset of diffraction (or, 
equivalently, the excitation of several photonic modes at 
the same frequency) is not restricted to the vicinity of 
the W or K/U points of the Brillouin zone and should 
be viewed as a more general characteristic of the optical 
response of the investigated photonic crystals. 

We notice that weak spectral features with a negative 
dispersion as a function of the incidence angle may also 
be recognized in the range 30° — 60° above the band gap 
peak (Fig. 6b and Fig. 6d). We do not pursue a quan- 
titative analysis of these features, whose identification 
is not straightforward. In fact photonic modes occuring 
at relatively high energy are more sensitive to disorder, 
which may broaden and mix fine spectral features, and 
also to dispersion of the refractive index. Further work 
is required in order to establish conclusively the presence 
of Van Hove singularities at energies above the band gap. 

IV. CONCLUSIONS 

In this paper a detailed study of the photonic band 
structure for artificial opals with close packed fee struc- 
ture is reported. The reduced photonic bands and density 



of states corresponding to a given value of the incidence 
angle in air are calculated in order to provide a guide 
for the interpretation of optical spectra. Two main fea- 
tures with opposite dispersions have been identified. A 
minimum of the reduced DOS corresponds to the stop 
band in the propagation direction, while maxima and van 
Hove singularities in the reduced DOS are associated with 
the excitation of photonic modes in the crystal due to 
diffraction by higher order planes in the photonic lattice. 
These modes depend on the investigated orientations of 
the sample with respect to the normal direction, reflect- 
ing the microscopic symmetries of the photonic lattice. 
By comparing with calculated reflection and diffraction 
spectra, it has been shown that van Hove singularities 
do indeed correspond to diffraction peaks and become 
also visible in reflectance spectra at large values of the 
incidence angle. 

The theoretical data are compared with microre- 
flectance spectra recorded at different incidence angles 
on bulk artificial opals. There is a very good agreement 
between theoretical and experimental data for what con- 
cerns the dispersion of the stop band. The predicted dis- 
persion of the excited photonic modes is validated by fea- 
tures observed in reflectance spectra for large incidence 
angles and different sample alignments along the symme- 
try directions LW and LK, thereby providing evidence 
for diffraction by higher-order planes (namely {111} and 
{200} families), or excitation of photonic modes out of 
the direction of light propagation. Finally, all the spec- 
tral features here discussed and measured in opals grown 
with different sphere diameter scale with the lattice con- 
stant, thus confirming their assignment to intrinsic fea- 
tures of the photonic band structure. The role of light 
polarization on the optical properties of opaline photonic 
crystals is currently under investigation. 
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